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ABSTRACT 

I examine the effectiveness of Kozai oscillations in the centres of galaxies and in par- 
ticular the Galactic centre using standard techniques from celestial mechanics. In 
particular, I study the effects of a stellar bulge pot ential and general rela tivity on 
Kozai oscillations, which are induced by stellar discs. iLockmann et al.l (|2008( ) recently 
suggested that Kozai oscillations induced by the two young massive stellar discs in 
the Galactic centre drives the orbits of the young stars to large eccentricity (e w 1). If 
some of these young eccentric stars are in binaries, they would be disrupted near peri- 
centre, leaving one star in a tight orbit around the central SMBH and producing the 
S-star population. I find that the spherical stellar bulge suppresses Kozai oscillations, 
when its enclosed mass inside of a test body is of order the mass in the stellar disc(s). 
Since the stellar bulge in the Galactic centre is much larger than the stellar discs, 
Kozai oscillations due to the stellar discs are likely suppressed. Whether Kozai oscilla- 
tions are induced from other nonspherical components to the potential (for instance, 
a flattened stellar bulge) is yet to be determined. 

Key words: Galaxy: centre ~ celestial mechanics 



1 INTRODUCTION 

Within 0.04 pc from Sgr A*, a swarm of young stars, 
known as the S- stars(Ec kart fc Genzell Il997l ) or SO-stars 
l|Ghez et al.ll2005h . orbit the central supermassive black hole 
(SMBH) in highly eccentric Keplerian orbits with random 
inclination, i.e. an isotropic distribution. These S-stars are 
typically of type B2 with a mass of ~ IOMq and a main 
sequence lifetime ~ 20 Myrs ((Alexanden |2005.) . Outside of 
these S-stars (0.04 pc to 0.4 pc), another population of 
young, massive stars, which consist of ~ 6 ± 1 Myr O- 
and early B-type stars, have a very different distribution. 
Rat her than being isotropic, these stars arc organized into 
one (|Lu et al.i .2006) or two disc(s) ( |Paumard et al. 200 6). 

The exact numb er of discs remains controversial. 
iPaumard etHI l|2006l ) find that these young stars are 
arranged in two stellar discs, whic h are inclined at ~ 
115 degr ees relative to one another (Paumard et al.| [200^ : 
iLevin fc Bc loborodov 2003,: iGenzel et al. 20031. The stars in 
one disc (which is the more massive disc) have roughly circu- 
lar orbits with the typical eccentricity < 0.4, while the stars 
in the other less massive disc are very eccentric (e > 0.6). On 
the other hand lliu et al.l l|2006l) find only one disc of stars, 
the more massive disc which IPaumard et al.l l|2006l ) iden- 
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tifies, and a "halo" population of young stars in random, 
highly inclined, and eccentric orbits. 

The origin of the young stars in th e central parsec o f the 
Galaxy remains an unsolved problem (|Alexandeill2005h . For 
the young stars between 0.04 and 0.4 pc, their arrangement 
in a disc suggests that they m ay have formed in-situ from 
the condensation of a gas d i sc (iLevin fc Beloborodov 20031: 



Navakshin fc Cuadral l2005l : iNavakshin. Cuadra fc Springell 
20071 : lLevinll2007l ). Such star formation would be expected 



in accretion discs at large radii due to the fragmentation 
from self-gravitv (|Paczvnskilll978l: iKolvchalov fc SvunvaevI 



trom selt-gravity llraczvnskilll^Yal: l-tvolvcnalov &: bvunvaevi 
1980l:lGoodinanll2003l : lThompson. Quataert fc Murravll2005l : 



Chansll200i r The more massive, more securely identified 
disc seem to favor this in-situ mode l as the stellar orbits 
are circular a nd dynamically cold (|Paumard et al.l l2006l : 
ILu et al]|2006l ). Alternatively, these stars may have formed 
in a large cluster at large radii, whi c h spiraled inward due to 
dynamical friction llGerhar dll200ll : E ansen fc Milosavlievid 
l2003l : iBerukoff fc Hansenll2006l ') . 

The problem of the S-stars' formation is even more 
daunting than that of the youn g, massive stars between 
0.04 and 0.4 pc (for a review see l Alexandeij|2005l ). For the 
S-stars, the tidal field from the SMBH and their current 
random orientation argues against both an in-situ disc for - 
mation scenario and a sinking cluster l|Ghez et al.l |2005| ) . 
The S stars may have for med via a diffe rent channel than 
the young, massive stars. iLevinI (|2003) has argued that 
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type 1/2 migration to small radii (iGoldreich fc Tremaine 
1978 ) , followed by resonant relaxation (|Rauch fc Tremaine 
1996 1 is a viable formation channel. On the other hand 



Perets. Hopman fc Alexanderl |20P7.) argues that the large 



population of massive perturbers near the Galactic centre 
(GC) increases the relaxation rate to such an extent that 
these S-stars may have been formed when binaries, which 
are scatter ed in to large ec centricities, are d isrupted near 
pericentre (|Hills|[T98& , Gould fc Quillenll2003D . 



Recently, a rather elegant scenario ha s been proposed 
by (|Lockmann. Baumgardt fc Kroupall2008l . hereafter LBK) 
for the origin of the S stars and their link to the young, 
massive stars which surround them. LBK showed that t he S- 
stars could be formed from ti dal disruption of binaries (|Hillsl 
1 19881 : iGould fc Quillenll2003l ) whi ch are driven to large ec - 
centricity from Kozai oscillations l|Kozailll962l : lLidovlll962l ) 
induced from the outer two young massive stellar discs. Ear- 
lier work by Subr, Karas and collaborators (|Subr et al.ll2004| ; 
ISubr fc Kar"3 l2005l : iKaras fc Subil |2007| ) studied the case 
for a single disc (with gasdynamic dissipative effects for 
the case of a fossil gas disc). This is indeed a very attrac- 
tive proposal as it naturally would explain their isotropic 
distribution, large eccentricity, and apparent youth. Using 
an N-body calculation wit h relativistic correction of up t o 
2.5 post-Newtonian orders (jLockmann fc Baumgardtll2008h . 
they showed that stars (or binaries) from the two discs can 
achieve eccentricities as large as e ~ 0.999. Relativistic cor- 
rections, which typicall y, can damp large eccentricities in 
the Kozai mechanism (IHolman. Touma fc Tremaind [l997l : 
iFabrvckv fc Tremaind |2007| ) appear not to be important 
(LBK). Binaries in these eccentric orbits that are disrupted 
near pericentre will result in one of the stars in a tight orbit 
around the SMBH. The reduction of semimajor axis from 
such a binary fission ev ent would be of order a factor of ten 
l|Gould fc Quillenll2003l ) and would be a viable mechanism 
for the formation of the S-stars. 



The beauty of LBK's scenario is that it ties in observed 
properties of the S-star distribution, namely their isotropy, 
eccentricity, and small semimajor axis with the observed 
properties of the young stellar discs, namely their high incli- 
nation of 115 degrees relative to each other. However, LBK's 
result is surprising as the Kozai mechanism is fairly delicate 
and can be greatly suppressed if additional perturbations 
to the gravitational potential such as a stellar cusp are in- 
cluded. Thus, I am motivated to examine of the basic physics 
of this scenario and ask the question: to what degree might 
Kozai oscillations be important in the centres of galaxies. 
In this paper, I study the nature of the Kozai mechanism 
central this scenario. I present my basic model and the basic 
equations for the Keplerian orbital parameters in ^ In SjS] I 
apply my basic model to the Galactic Centre. I numerically 
compute the evolution of a star subject to the perturbed po- 
tential from a single disc, a single disc with a stellar bulge, 
two discs, and two discs with a bulge. I show that including 
the spherical stellar distribution of sufficient mass suppresses 
the Kozai mechanism. I then argue that the spherical stellar 
distribution in Galactic centre is more than adequate of this 
purpose. I discuss some of these implications and conclude 
in 



2 SECULAR EVOLUTION 

I now describe the model problem central to this study. Two 
massive stellar discs with mass Mi and A'h orbit a central 
SMBH with mass Mo 3> Mi , M2 at an inclination relative 
to one another of iq. The surface density of these two stellar 
discs. El, 2, is assumed to be a power law with index —p, 
i.e.. El, 2 ~ 'Si,2,o{r /ro)~^ , where r is the radial coordinate, 
ro is a reference radius, and Ei,2.o is the normalization for 
disc 1 and 2 respectively. Without loss of generality, I orient 
my axis such that the reference plane is in the plane of the 
Ml disc and I presume that Mi ^ M2. Also surrounding 
the SMBH is a spherical stellar power law distribution with 
index — g, n* — no{r /ro)~'' ■ A test body orbits the SMBH in 
a near-Keplerian orbit with semimajor axis, a, inclination, 
i, and eccentricity, e. 

I first consider the case of a perturbing mass, Sm, in a 
circular orbit of radius, r, around the central mass and its 
effect on the test body. I will work exclusively in the secular 
approximation. This pr oblem is well studied by many 
auth or s (see for in s tance llnnanen et al.|[l997l : Kiseleva et al 
1 19981: iFord et all |2000|'). In particular. lEggleton et al 
(1 19 981) (see also [Eggleton fc Kiseleva-EggletonI 12001 



Fa brvckv fc Tremainel 120071 ') has introduced a nice formal- 
ism which I find very flexible, powerful, and especially 
useful for studying multiple perturbing bodies which 
are at large in clinations relative to one ano t her. Using 
the notation o f Egglcton fc Kiseleva-EggletonI (|200lh and 
IFabrvckv fc Tr emainc (2007) and dropping dissipative and 
tidal terms, the governing equations are: 



1 de 
e dt 



- (1 - e") ~ (4S,e - Sgg) q + Sqhh] , (1) 

~ = (1 - e') S,he - (4e' + l) SehQ + de^Segh, (2) 

where e is the Laplace- Runge-Lenz vector, whose magni- 
tude is the eccentricity, e, h is the reduced orbital angular 
momentum vector, h and e are the normalized vectors in 
the direction of h and e respectively, and q — h x e is the 
normal vector that completes the triad. {h,e,q). The dis- 
turbing tensoiQ is Sxy ~ C \5xy — 3{h ■ x){h ■ y)j , where 

h is the normalized angular momentum vectors of the per- 
turbing mass and the constant, C, is 



C = ■ 



-1/2 



(3) 



where n{a) — \J GMo /a^ is the mean motion. Finally Zgk 
and represents the apsidal motion from the effects of 
general relativity and the ste llar cusp. The precession term 
due to general relativitv is ([ E ggleton fc Kiseleva-EggletonI 



I2OOII : IFabrvckv fc Tremaindl2007l ) 



7 3 rg 



2 a 1 



(4) 



where r-g = 2GMq /(? is the gravitational radius. The preces- 
sion due to the stellar bulge potential is Ijlvanov et al.ll2005l . 
see their eq.(14) and (15) and appendix A) 



^ I am unaware of a name for this tensor, so for the sake of 
nomenclature, I have chosen to call this the disturbing tensor. 
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M«(a) 



(5) 



where Mt(a) = 47r mtn{r)r^ dr is the mass of the bulge 
stars inside of sphere whose radius is equal to the test body's 
semimajor axis, a, 

r(5/2 - q) 

is a constant of order unity, and F is the gamma function. 

I now calculate the effect of a disc. Taking 5m — > dm = 
2n'Erdr, 1 find that the effect of a disc is a modification of 
equation ([HJ to be 



Cd = n(a)Mo a 



-1/2 



2-K 



drT,r dr. 



(7) 



where Cd is the constant associated with the disc, rin is the 
disc inner radius, and rout is the outer radius of the disc. 
For rout > rin, this gives 

3 



Cd = n{a) 



off 



Mo 



1 



-1/2 



(8) 



where M^b ~ 2/(f + p)7rEi,oro is the effective mass of the 
disc, and Md — 2n J'^""' drrE is the mass of the disc. Note 
that the effect of a massive extended stellar disc with rout 
rin can be reduced to a single ring with mass Moff at a 
radius of ri,2in for p > — f. For p — — f, the contribution 
is logarithmic in radius and for p < — f, the outer radius is 
more relevant. 

The use of the governing equations ((T| and Q is espe- 
cially advantageous for studying the secular effects of mul- 
tiple pertu rbing rings. Like the Lagran g e-Laplace planetary 
equa tions (iBrouwer fc Clemencd Il96ll : iMurrav fc DermottI 
I2OOOI ). the governing equations ^ and ([2]) are linear. Hence, 
the effect of additional perturbing masses is just a matter of 
adding their associated disturbing function (in the case of 
the Lagrange-Laplace planetary equations) or the perturb- 
ing potential, Sxy, in the case of equations ([TJ and ((2}. The 
calculation of the inclined disturbing function in the case 
of the sec ularly averaged Laplace-L a ngrange's planetar 



ot tne sec ularly averaged l^aplace-L a ngrange s planetary 
equations (iBrouwer Sz Clemencel 196ll: Innanen et al.l 19971: 



Kiseleva et al.lll998l : iMurrav fc DermottI I2OO0I : I Ivanov et al.l 
20051 ) is rather involved and I have found it to be numerically 
difflcult to solve. By contrast including additional masses is 
trivial for the governing equations ((T| and 0. For a number 
of perturbing masses, 5mi, the disturbing tensor, Sxy, is a 
sum over all disturbing tensors from all the masses, or 



(9) 



where hi is the normalized angular momentum vectors of 
the perturbing mass Smi. The associated constant, d, is 

3 



, ,Smi ( a 
Ci = n(a)—— — 



1 



-1/2 



(10) 



where ri is the radial position of the perturbing mass. 



APPLICATION TO THE GALACTIC 
CENTRE 



In the case of the Galactic centre, (|Paumard et al.l |2006| ) 
find that the stellar distribution follows p ~ 2, ri,2,in = 0.1 




Figure 1. Simple illustration of the Kozai oscillation for one disc 
of mass, A/off = 2 X IO^Mq for a test particle (star) with an 
initial inclination of 60 degrees relative to the plane of the disc. 
The test particle's eccentricity and inclination are shown as a 
solid and dotted lines respectively. The precession from general 
relativity and the stellar bulge potential are ignored in this case. 



pc, and ri,2,out = 0.5 pc for the two disc, which are oriented 
at 115 ± 7 degrees with respect to one another. The upper 
limit for the mass of the clockwise disc is Md,i < IO^Mq, 
and the upper limit for th e mass of the counte rclockwise 
disc is Md,2 < 5 X lO^Mp jP aumard et al.ir2006l ). which is 
consistent with lNavakshin et al . (2006) 's dynamical limits of 
< IO^Mq for the mas ses of two discs. For a t op-heavy initial 
mass function (IMFl. lPaumard et all l|2006l ) derives a lower 
limit for the stellar mass of the two discs to be 35OOM0 and 
14OOM0. Hence for a range in disc masses between Md,i — 
3500 - 1O*M0 and Md,2 = 1400 - 5OOOM0, the effective 
disc masses are Mi, off ~ 0.6 — 2 x IO^Mq and M2,cff ~ 
0.3 - 1 X IO^Mq. 

1 rmmerically calculate the evolution of equa tion ^ and 
I2I), using a standard Runge-Kutta algorithm (jPress et al.l 
1 1992 ). As an illustration, I plot the eccentricity (solid line) 
and inclination (dotted line) for a test particle that is in 
orbit around a SMBH and is initially inclined relative at 
60 degrees to a single massive disc with M^g = 2OOOM0 in 
Figure [1] For this case, I have taken Zen ~ and — 
0. Note the characteristic behavior of the Kozai oscillation, 
where the inclination and eccentricity vary in phase and the 
initial inclination determines the amplitude of the oscillation 
due to the conservation of the Kozai integral or z-com ponent 



of th e angular momentum vector, Lz — (1 — e^) cos^ i (|Kozail 
Il962h . 

The axissymmetry, but not spherical symmetry, of the 
potential conserves Lz , but not the total angular momentum 
vector, L. Thus, an initially inclined low eccentricity orbit 
achieves very hig h eccentricities. For th e cas e of the region 
aroun d a SMBH, ISubr fc Karag (|2005l ) and iKaras fc gubil 
l|2007l ) studied the effect of fossil gas disc on the orbits 
of a nuclear star cluster. In the case of the a single disc. 
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ISubr fc Kar"a3 (|2005l ') and lKaras fc Subij l|2007l ) showed that 
there are two kinds of orbits: orbits which librate around 
uj = n/2 and 3tt/2 and orbits which span over uj = [0, 2n], 
where co is the argument of pericentre. In the particular 
case of a fossil gas disc, dissipative interactions from star- 
crossings of the fossil gas disc dissipa tes energy, resulting in 
a slow decay of the semi-major axis l|Karas fc SubjlioOTl '). 

I now study the effect of Kozai oscillations including 
the effects of general relativity, Zgr and the second disc. In 
agreement with the LBK's claim, I do not find the effect of 
relativity to be significant for eccentricities up to 0.999. 

I now include the effect of the second disc, which is 
less massive at « 1.4 — 5 x IO'^Mq than the first disc 
(3.5 — 10 X IO^Mq). The observed inclination is fairly nar- 
row 115 ± 7 degrees, so I choose to fix the inclination at 
115 degrees. In Figure [21 I show the maximum eccentricity, 
fimax, reached as a function of initial cos i relative to the ref- 
erence plane (which is the plane of the disc for the one disc 
case and the plane of the more massive disc in the two disc 
case). For two discs, Cmax as a function of cosi is consider- 
ably more complex than the case for one disc. Note that for 
low mass discs (either for one disc or two disc case), there 
is insufficient time (7 Myrs) for stars to reach high eccen- 
tricity. For two discs, the range of initial inclinations that 
generate large eccentricities is larger compared to the sin- 
gle disc. Also note, that large eccentricities can be reached 
for two discs in the neighborhood around the inclination of 
the second disc. The addition of a second disc increases the 
available inclinations over which test bodies, i.e., stars, can 
reach large eccentricities. Note also that even in the most 
optimistic scenario, stars do not oscillate to high eccentric- 
ity for cosi ~ ±1. Stars that are close the to more massive 
stellar disc do not reach large inclination in spite of the sec- 
ond disc. Hence, stars and binaries that are driven to large 
eccentricities must come from the less massive disc, which 
supports LBK's result that stars in the CCW disc (the less 
massive disc) are more likely to be driven to large eccentric- 
ity. 

The axissymmetry of a single disc conserv es L^. As a re- 
sult, a polar plot of e-ui (see for instanc dSubr fc KaradlioOSl : 
iKaras fc SubillioOTi ) of an orbit generates a closed curve. On 
the other hand, when a second disc is included, axissymme- 
try is broken and an orbit no longer generates close curves 
in a polar e-ui plot. 

When the stellar bulge potential is included, Z*, the 
dynamics changes completely. For the bulge potential I take 
g = 1.4 (which gives k « 0.6 for eq.[5]) and no = po/M© 
is the nu mber density of stars, where po is the mass density 
of stars (jGenzel et al.ll2003l : IYu et al.ll2007l 'l. I plot its effect 
as a function of M, (a — 0.1 pc) in Figure [3] for the most 
optimistic scenario Md,i = lO^M© and Md,2 = 5OOOM0. As 
this figure shows, sufficiently large stellar bulges (M*(a = 
0.1 pc) > 35OOM0 for one disc and (M.(o = 0.1 pc) > 
45OOM0 for two discs) suppresses the Kozai mechanism. 

The mass found that is needed to suppress Kozai oscilla- 
tions agrees well with simple analytic arguments. The Kozai 
mechanism is known to be suppressed by non-Keplerian 
contributions to the potential such as those introduced by 
general relativity or addition al masses in the system (se e 
the case for multiple planets, iTremaine fc Zakamskall2004h . 
Typically if the period for apsidal precession is of order or 
shorter than the period of the Kozai oscillations, then the 




cos(i) 

Figure 2. Maximum eccentricity as a function of initial inclina- 
tion for one disc (dashed lines) with Mi off = 0.6 (lower dashed 
line), 1.2 (middle thick dashed line), 2 X 10^ Mq (upper dashed 
line) and two discs with Mi eff = 0.6 (lower solid line), 1.2 (mid- 
dle thick solid line), 2 X IO^Mq (upper solid line) with respective 
M2 cff = 0.3, 0.6, 1 X IO'^Mq). As a function of phase space, note 
that the regions available for large eccentricities (e > 0.95) is 
much larger for two discs as opposed to one disc. The vertical 
dotted line indicates the position of the second disc relative to 
the first for the parameters of the two discs in the Galactic cen- 
tre, i.e., an inclination of 115 degrees. 

Kozai oscillations are suppresse d (jTremaine fc Zakamskal 
I2OO4I : iFabrvckv fc TremainellioOTi ) . The period of the Kozai 
oscillation is of order (from studying eq.[3]): 



ri,2.i 



Mo 



(11) 



whereas the timescale for apsidal precession due to the stel 
lar bulge is of order (from studying eq.[5]) 

-1 M.(a) 



Mo 



(12) 



Hence, the requirement that the apsidal precession period 
be longer than the period for Kozai oscillations imply M, < 
Moff , matching the expectations from the more detailed cal- 
culation in Figure [3l 

I now argue that the mass of stellar bulge is much 
larger than what is needed to suppress the Kozai mecha- 
nism. From observations, the measured M«(a = 0.1 pc) is 
~ 6 X 10 *M(7), using the valu es for the central stellar den- 
sity from lGenzel et al.l (|2QQ3| ) (see also Schodel et al. 2007^). 
The observed mass is over an order of magnitude larger than 
what is needed to suppress the Kozai mechanism. It is sig- 
nificantly larger than the mass of the two young stellar discs 
and, therefore, Kozai oscillations induced by the two young 
massive stellar discs are likely suppressed. 

The suppression of Kozai oscillations due to a spherical 
distri bution of stars is well known for the c a se of a single 
disc (llvanov et al.ll2005l : iKaras fc Subil 120071 ). Ilvanov et all 
(I2OO5I 7 argues that in the case of a single ring, Kozai oscilla- 
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oretical grounds, such a significantly fiattened bulge may 
be un likely because resonant relaxation ( Ranch fc Tremaing 
Il996l ) would isotropise the bulge stars ( Levinll2007l . see for 
instance,) over their 10 Gyr lifetime. Observationally, there 
is no evidence for a flatten bulge as the velocity distribution 
of the late-typ e stars appear to be consistent with isotropy 
l|Genzel et al.| [r996V 

The Kozai mechanism may be more applicable in other 
galactic nuclei such as M31, which has a mass ive stellar 
disk (the P1/P2 disk). In M31, the P1/P2 disc (|Tremaind 
ll995l : |Peiris fc Tremainell2003l : IChang et al.ll2007l) is ~ 10% 
of th e mass of the 1.4 x 10* M© SMBH l|Bender et all 
l2005f ). whereas the stellar bulge potenti al is much smaller, 
i.e., M « < lO'^Mp) at 1" or « 4 pc jPeiris fc Tremaind 
l2003f ). IChang et al] (|2007l ') has suggested that the non- 
axisymmetric potential of the P1/P2 disc modifies gas or- 
bits such that they are confined to be inside of 1 pc around 
the SMBH. High inclination stellar orbits may also undergo 
Kozai oscillations in this case. The implications for these 
stars undergoing Kozai oscillations in the nucleus of M31 
would be an interesting topic for further study. 



Figure 3. Maximum eccentricity as a function of stellar mass 
enclosed for a = 0.1 pc for one disc (dashed line, Mi efj = 
2 X lO^M©) and two discs (solid line, Mi j.ff = 2 X lO^M© and 
M2 eff = IO^Mq). Also shown is Mi eff and M2,eS for the com- 
parison of scales. For sufficiently large stellar bulges (Mt(a = 
0.1 pc) > 3500Mq for one disc and M,{a = 0.1 pc) > 4500Mq 
for two discs), large eccentricities cannot be achieved and the 
Kozai mechanism is suppressed. 



tions are suppressed for initially low eccentricity orbits for a 
sufficiently massive spherical stellar distribution. Similarly, 
iKaras fc §ubil l|2007l ) also showed that a spherical stellar dis- 
tribution whose mass is comparable to the perturbing disc 
mass will suppress Kozai oscillations for initially circular or- 
bits (see the third panel of their Figure 2). The results of 
this paper are in broad agreement with these previous re- 
sults and point out that inclusion of additional rings, which 
breaks axial symmetry, do not change the basic results for 
small initial eccentricities. 



4 DISCUSSION AND CONCLUSION 

I find that the apsidal precession induced by the stellar bulge 
(using a realistic estimate for its mass) greatly reduces the 
impact of the Kozai mechanism in the GC. Hence, the Kozai 
oscillations central to LBK's elegant mechanism is likely not 
due to the stellar discs. The calculation in LBK is more re- 
alistic in that it captures the dynamics of the system with- 
out resorting to the perturbative scheme used in the present 
study. However, this study is complementary because it out- 
lines the regions of parameter space where secular effects are 
important. 

According to the present work, the young stellar disc is 
not likely to induce Kozai oscillations, but other nonspher- 
ical components to the potential may able to do so if they 
are sufficiently strong. One possibility is a significantly flat- 
ten (of order unity) stellar bulge. A detailed study of the 
degree of flattening required to induce Kozai oscillation is 
interesting, but it is beyond the scope of this work. On the- 
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